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0.  Abstract  (Continued) 


)  minimize  the  power  in  the  specified  sidelobe  region  subject  to  the  con¬ 
straint  that  the  array  response  in  the  mainbeam  direction  is  held  constant. 
It  is  shown  that  forms  2  and  3  of  the  minimization  problem  lead  to  identical 
results,  and  that  the  results  obtained  with  form  1  are  numerically  close  to 
those  obtained  with  forms  2  and  3.  We  conclude  that  if  all  except  one  de¬ 
grees  of  freedom  of  the  element  weights  are  utilized  to  minimize  pattern 
power  in  a  specified  sidelobe  sector,  with  the  remaining  one  degree  of 
freedom  used  to  satisfy  a  mainbeam  constraint  condition,  then  a  pattern 
can  be  formed  with  low  average  sidelobe  power  in  a  broad  sector.  The 
larger  the  number  of  array  elements,  the  lower  will  be  the  average  power 
in  the  specified  sidelobe  sector.  The  price  paid,  however,  for  this  sector 
sidelobe  power  reduction  is  a  marked  degradation  of  the  remainder  of  die 
pattern. 
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1.  INTRODUCTION 

The  growing  interest  in  adaptive  antennas  makes  it  worthwhile  to  examine  the 
theoretical  limits  to  what  can  be  achieved  with  reduction  of  sidelobe  levels  of  array 
antenna  patterns.  In  tills  report  we  address  the  problem  of  determining  the  mini¬ 
mum  possible  power  that  can  be  present  in  a  specified  sidelobe  region  subject  to 
certain  constraints.  We  consider  three  different  forms  of  the  minimization  problem: 

1)  minimize  the  ratio  of  the  power  in  the  specified  sidelobe  region  to  the  power 
in  the  main  beam  region: 

2)  minimize  the  ratio  of  tile  power  in  the  specified  sidelobe  region  to  the  power 
in  the  main  beam  direction;  and 

3)  minimize  the  power  in  the  specified  sidelobe  region  subject  to  the  constraint 
that  the  array  response  in  the  main  beam  direction  is  held  constant. 

When  we  refer  to  the  main  beam  region  or  main  beam  direction  we  are  refer¬ 
ring  to  a  pattern  prior  to  the  indicated  minimization  since  the  main  beam,  in  gen¬ 
eral,  will  shift  or  transform  in  shape  as  a  result  of  the  change  in  array  weights. 

It  is  shown  that  forms  2  and  3  of  the  minimization  problem  lead  to  identical 
results,  and  that  the  results  obtained  with  form  1  are  numerically  close  to  those 
obtained  with  forma  2  and  3. 

(Received  for  publication  12  November  1981) 


For  all  three  forms  of  the  sidelobe  power  minimization  problem  we  consider 
an  equispaced  linear  antenna  array  consisting  of  2N  +  1  isotropic  elements  spaced 
a  distance  d  from  each  other  (see  Figure  1).  The  field  radiation  pattern  of  the 
array  is 

N 

F(u)  =  I  a^™ 
n=-N 

where 

an  «  the  n^  complex  array  element  weight 


with 


k 


2  T 
\ 


and 


d  =  the  pattern  angle  measured  from  broadside. 
The  corresponding  power  pattern  is  then  given  by 

N  N 

F{u)  2  *  T  Z  a  a  e'j(n'm)u  . 

w  m  n 

n=-N  n- -N 


(1) 


l'he  direction  of  the  ma.n  beam  prior  to  adjusting  the  array  coefficients  to  reduce 
tnc  power  in  a  specified  side-lobe  region  will  be  assumed  throughout  to  be  the 
broadside  direction,  u  =  d,  but  the  analysis  presented  here  can  be  readily  extended 
o  cover  any  specified  main  beam  direction. 

2.1  Minimize  tin-  Ratio  of  file  Power  in  the  Specified  Sidelobe 
Region  to  (lie  Power  in  the  Main  Beam  Region 

The  array  pattern  corresponding  to  uniform  amplitude  and  zero  phase  element 
coefficients  is 


.  f 2N  +  1  \ 

1  stn(-2— )u 

2N  +  1  .  u 

Sln  2 

We  take  the  main  beam  region  to  be  half  the  interval  between  the  first  nulls  on 
either  side  of  broadside,  or 

"  UB  -  uS  UB 
where 


u 


_  -  £ 

B  ‘  2N  +  ! 
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The  power  in  the  main  beam  region,  Pg,  is  then  given  by 


PB  =  iff  /  lF(u)!2  du  . 


‘B 


Substituting  Eq.  (1)  for  |F(u)|  and  performing  the  integration  then  gives 


N  N 

’b  =  kd  £  £ 

m  =  -N  n  =  -N 


a’  a 
m  n 


sin  [  (n-m)  Ug] 
(n-m) 


2  u 


N  N 


=  ~kT  ^  ^  am  an  sinc  1  (n'm)  UB]  ‘ 


(2) 


m=-N  n=-N 


The  sidelobe  region  is  specified  by 

u  -  £  ^  u  £  u  +g 
o  o 

where 

uQ  *  kd  sin  6Q  . 

The  power  in  the  sidelobe  region  is  then  given  by 
uo+e 


ps  •  A  /  ipfa>i2<ta 

uQ-e 

N  N  “0+e 

•H  2  2  *m*n  / 

m*-N  n*-N  u  -e 

o 

2e  ^  «  j(n-m)  u 

■  M  2  I  am  an  e  0  sine  [  <n-m)«J  . 

m*-N  n*-N 


(3) 
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We  now  desire  to  minimize  the  ratio  of  the  power  Pg  in  the  Bidelobe  region  to 
the  power  Pfi  in  the  main  beam  region.  Thus  we  must  minimize  the  quantity 
M  =  PS/PB-  or  substituting  from  Eqs.  (2)  and  (3). 

N  N 

I  S  am  an  Hmn 

M  *  r  — N~  -  (4) 

£  ^  am  an  Gmn 

m*-N  n»-N 

where 


r 


(5) 


e 

UB  2N+1 

~j  (n-m)  u 

Hmn  =  *  °  sine  f  (n-m)  t] 

and 


Gmn  ’  ainc  1  <n'm)  UB>  ’ 

Equation  (4)  for  u  can  be  written  in  the  form 

a+  H  a 
M< a)  ■  r  r — ~ 
a  G  a 

where  a  is  the  vector  of  complex  array  element  coefficients 
a-(N- 1)  I 


aN-l 
.  aN 


and 


+ 


a 


is  the  Hermitian  transpose  of  a. 


(6) 


(7) 
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For  p  to  be  a  minimum  it  is  necessary  for  m  to  be  stationary  with  respect  to 
small  variations  in  a;  that  is,  p(a  +  6a)  =  p(a)  to  the  first  order  for  any  small 
variation  6_a.  Equating  p( a  +  6a)  =  p(a)  in  Eq,  (7)  leads  to  the  eigenvalue 
equation* 

Ha  =  pGa 


or 


Ha  =  p'  Ga,  (8) 

letting  n'  =  n/r. 

Now  both  H  and  G  are  Hermitian  matrices  and  G  is  positive  definite  since 

a  Ga  equals  the  power  in  a  region  of  the  array  pattern  and  so,  is  greater  than  zero 

2 

for  any  set  of  array  coefficients  not  identically  equal  to  zero.  Hence  the  eigenvalue 
equation  Eq.  (8)  has  2N  +  1  linearly  independent  eigenvectors  a;  i  =  1,  2, .  .  .  ,  2N  +  1, 
satisfying  the  orthogonality  relation 

a!  Ga  j  =  1,  i  =  j 

=0,  i  t  i  . 

The  smallest  eigenvalue  of  Eq.  (8)  when  multiplied  by  r  as  defined  by  Eq.  (5)  is 
the  desired  minimum  value  of  the  ratio  of  the  power  in  the  sidelobe  region  to  the 
power  in  the  main  beam  region. 

2.2  Minimize  the  Ratio  of  the  Power  in  a  Specified  Sidelobe 
Region  to  the  Power  Density  in  the  Direction  of  the 
Peak  of  the  Main  Beam 

This  minimization  is  obviously  equivalent  to  maximizing  the  ratio  of  the 
power  density  in  the  main  beam  direction  (that  is,  broadside)  to  the  power  in  the 
specified  sidelobe  region.  Now  the  power  density  in  the  direction  u  =  0  is  given 

by 


1.  Pease,  M.  C.  ,  III  (1985)  Methods  of  Linear  Algebra,  Academic  Press,  New 

York,  p.  91. 

2.  Gantmacher,  F.  (1960)  The  Theory  of  Matrices,  Vol.  1,  Chelsea  Publishing 

Co. ,  New  York,  p.  338. 
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N  N 

£  E  am  an 
m»-N  n“-N 


+  . 

a  A  a 


where  a  is  the  vector  of  element  coefficients  and 


11.. 


HI 

K  '  ' 

!  *'i 


As  in  the  analysis  of  the  first  form  of  the  minimization  problem,  the  power  in  the 
side-lobe  region  is  given  by 

N  N 

ps  -  sir  I  l  *S,  »m„  »„ 

m  =  -N  n=-N 


=  kl  a.+%  • 


with  the  matrix  H  defined  by  Eq.  (6).  Hence  we  wish  to  maximize  the  ratio 


a  A  a 

2e,  +  „ 
kd  -  H  - 


As  above  we  are  again  led  to  an  eigenvalue  problem;  namely. 


=  td  "  H  5. 


In  this  case,  however,  because  of  the  fact  that  A  is  given  by  Eq.  (9)  there  is  only 
one  non -zero  eigenvalue3  given  by 


L  l  J 


with  a  corresponding  eigenvector 


Ll  J 


(12) 


Note  the  n  equals  the  sum  of  all  the  elements  of  H  1  divided  by  -jA  and  that  the 

- 1  Ka  2  e 

array  coefficients  are  the  sums  of  the  respective  rows  of  H  divided  by 


2.3  Minimize  the  Power  in  the  Specified  Sidelobe  Region  Subject 
to  the  Conitraint  that  the  Value  of  the  Array  Pattern  in  the 
Direction  of  the  Peak  of  the  Main  Beam  it  Held  Constant 

The  field  radiation  pattern  of  the  array  in  the  direction  u  =  0  is  given  by 


F(0) 


N 

I  * 

n*-N 


so  that  the  constraint  that  F(0)  should  remain  constant  is  equivalent  to  requiring 
that 


c+  a  *  F(0) 


constant 


(13) 


3.  Cheng,  D.K.  and  Tseng,  F.  I,  (1965)  Gain  optimization  for  arbitrary  antenna 
arrays,  IEEE  Trans.  Ant.  Prop. ,  Vol.  AP-13,  pp.  973-974. 


A 


12 


where 


1 

1 

c  * 

.  1 

Since  the  power  in  the  sidelobe  region  is  given  by  Eq.  (10), 

P  *  a+  H  a 
S  kd  -  - 


the  problem  here  is  to  determine  the  vector  a  of  array  element  coefficients  that 
minimizes  a+Ha  subject  to  the  constraint  Eq.  (13).  The  solution  of  the  constrained 

—  —  4 

minimization  problem  is  given  by 


H  1  c 

a  *  — - f-  F(0)  .  (H) 

c  H  c 

The  denominator  of  this  expression  is  the  sum  of  the  elements  of  H  1  and  the  i**1 
element  of  a  is  the  sum  of  the  elements  of  the  i**1  row  of  H  1  divided  by  the  sum  of 
all  the  elements  of  H*  *.  Hence  the  sum  of  the  element  coefficients  is  F(0)  as 
desired. 

The  minimum  sidelobe  region  power  is  given  by 


u 

kd 

a+  H  a 

Is, 

(  h'*£F(o)1+  h 

[H^c  F(0» 

kd 

<£+ 

H’1 

c)2 

2 i 

|  F<0)|  2  c+ 

H"1 

c 

kd 

(c+  H”1 

c)2 

2e 

|f(o)J  2 

kd 

+„-i/ 

C  n  C 

4.  Frost,  O.  L. ,  III  (1970)  Adaptive  Least  Squares  Optimization  Subject  to  Linear 
Equality  Constraints.  SEL-70-055,  'tech.  Rep.  No.  6796-2,  Information 
Systems  Laboratory,  Stanford  University,  California,  p.  6. 
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where,  we  have  made  use  of  the  fact  that  H  and  hence  H  *  is  Hermitian  to  set 
(H"  1)+  =  H"1,  For  this  value  of  the  minimum  sidelobe  region  power,  the  ratio 
of  the  power  density  at  u  *  0  to  the  sidelobe  region  power  is 


Note  that  this  ratio  is  exactly  the  value  found  above,  Eq.  <  11).  in  the  analysis  of 
the  second  form  of  (he  minimization  problem,  for  the  maximum  ratio  U  of  the 
power  density  at  u  *  0  to  the  power  in  the  specified  sidelobe  region.  Hence  the 
minimum  sidelobe  region  power  obtainable  with  array  coefficients  that  satisfy  the 
constraint  Eq.  (13)  is  the  same  sidelobe  region  power  that  minimizes  the  ratio  of 
sidelobe  power  to  the  power  density  at  u  =  0. 

Comparing  the  expression  for  the  array  element  coefficient  vector  a  found  in 
Eq.  (14)  with  that  found  in  the  second  form  of  the  minimization  problem.  Eq.  (12), 
we  see  that  the  two  expressions  agree  to  within  a  multiplicative  constant.  This  is 
effectively  equality  since  eigenvectors  are  determined  only  to  within  an  arbitrary 
multiplicative  factor.  Thus  forms  2  and  3  of  the  minimization  problem  lead  to 
identical  solutions. 


I 


•Si 


3.  NUMERICAL  RESULTS  AND  DISCUSSION 

FORTRAN  computer  programs  were  written  and  run  on  a  CDC  6600  to  calcu¬ 
late  the  element  weights  and  corresponding  power  ratios  and  power  patterns  for 
forms  1  and  2  of  the  minimization  problem.  The  inputs  were  the  number  of  ele¬ 
ments  in  the  array,  2N  +  1,  the  element  spacing,  d,  in  terms  of  the  wavelength 
(that  is,  half  wavelength  spacing  is  given  by  d  *  1/2),  the  center  of  the  sidelobe 
region  in  degrees,  0o»  and  the  upper  half  width  of  the  sidelobe  region  in  degrees, 
£0.  In  terms  of  the  inputs  the  sidelobe  region  in  "u  space"  is  given  by 

u  -eSuSu+e  (15) 

o  o 

where 

u  =  2ird  sin  0 
o  o 
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and 


e  ■  2vd  (sin  (0Q  +  AO)  -  sin  flQl 

Note  that  the  sidelobe  region  defined  by  Eq.  (15)  is  symmetrical  with  respect  to 
die  center  of  the  region  in  "u  space"  but  not  in  "0  space". 

The  solution  of  the  generalized  eigenvalue  Eq.  (8)  involved  in  the  first  form 
of  the  minimization  problem  was  solved  using  the  IMSL  program  EIGZC  and  the 
inversion  of  the  matrix  H  in  Eq.  (12)  in  the  second  form  of  the  minimization  prob¬ 
lem  was  performed  using  the  IMSL  program  LEQ2C.  For  all  calculations  the 
sidelobe  center  0Q  taken  equal  to  35°  and  the  element  spading  d  was  set  equal 
to  l/2.  For  N  ■  3  (7  elements)  calculations  were  performed  for  AO  =  15°,  25°, 
35°,  45°.  and  55°;  and  for  N  =  4  and  5  (9  and  11  elements,  respectively)  calcula¬ 
tions  were  performed  for  AO  =  45°  and  55°. 

Table  1  gives  the  lower  and  upper  limits  of  the  sidelobe  region  in  degrees, 

0 _  and  0+,  for  the  various  values  of  AO  used. 


Table  1.  Lower  (e  )  and  Upper 
(0+)  Limits  of  the  Sidelobe  Re¬ 
gion  Corresponding  to  a  Sidelobe 
Region  Center  of  35°  and  an 
Upper  Half  Width  of  A© 


(o) 

AO 

(o) 

0_ 

(o) 

*+ 

15 

22.40 

50 

25 

16.33 

60 

35 

11.97 

70 

45 

9.34 

80 

55 

8.46 

90 

Table  2  gives  the  ratio  of  sidelobe  power  to  main  beam  region  power,  Pg/Pg* 
for  the  first  form  of  the  minimization  problem,  and  the  ratio  of  the  sidelobe  power 
to  the  power  in  the  broadside  direction,  Pg/|F(0)l^,  for  the  second  form  of  the 
minimization  problem.  For  comparison  purposes  the  same  quantitities  are  also 
tabulated  for  a  uniformly  weighted  array.  Note  that  as  die  number  of  elements  in¬ 
creases  the  sidelobe  power  decreases,  since  more  nulls  are  available  to  be  placed 
in  the  sidelobe  region. 


15 


Table  2.  Ratio  of  Sidelobe  Power  to  Main  Beam  Region  Power,  Pq/Pb.  for  the 
First  Form  of  the  Minimization  Problem,  and  Ratio  of  Sidelobe  Power  to  Broadside 
Direction  Power,  Ps/|F(0)|  ,  for  the  Second  Form  of  the  Minimization  Problem. 
Also  shown  for  comparison  are  the  same  quantities  for  a  uniform  array.  Numbers 
in  parentheses  are  exponents  of  10 


Seven  Elements 

(o) 

Ad 

15 

25 

35 

45 

55 

VPB 

0.402(-8) 

0.  1 17 (-5) 

0.  277  (-4) 

0.  151(-3) 

0.  259(-3) 

uniform 

0.  397(-l) 

0.  515(-1) 

0.  659(-l) 

0.  111(0) 

0.  136(0) 

ps/|f(o)|2 

0.  1 92 ( -8) 

0.  606(-6) 

0.  159(-4) 

0.  949(-4) 

0.  169(-3) 

uniform 

0.  881(-2) 

0.  114(-1) 

0.  146(- 1) 

0. 246(- 1) 

0.  301(-1) 

Nine  Elements 

(o) 

AO 

45 

55 

*VPB 

0.  795(-5) 

0.  164(-4) 

uniform 

0.  707(-l) 

0.  843{-l) 

ps/!f(o)|2 

0.  464(-5) 

0.  101(-4) 

uniform 

0.  12  2  (-  1) 

0.  145(- 1) 

Eleven  Elements 

(o) 

Ad 

45 

55 

fyPB 

0.  416(-6) 

0.  104(-5) 

uniform 

0.  607 (-1) 

0.  648(-l) 

ps/|f(0)|2 

0.  225(-6) 

0.597{-6) 

uniform 

0.  855(-2) 

0.  913(-2) 
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The  element  weights  for  the  first  and  second  forms  of  the  minimization  prob¬ 
lem  were  found  to  be  surprisingly  close  to  one  another,  differing  by,  at  most, 
about  four  percent  for  the  cases  studied.  The  power  patterns,  when  plotted  for 
corresponding  cases  of  the  two  forms  of  the  minimization  problem,  are  virtually 
indistinguishable  from  each  other.  In  Figure  2(a-e)  we  show  the  seven  element 
patterns  obtained  for  the  second  form  of  the  minimization  problem,  in  Figure  3(a 
and  b)  we  show  the  nine  element  patterns,  and  in  Figure  4{a  and  b)  the  eleven  ele¬ 
ment  patterns.  For  reference  purpose  the  uniform -weighted  seven,  nine,  and 
eleven  element  patterns  are  shown  in  Figures  5a,  5b,  and  5c.  All  patterns  are 
normalized  to  have  the  value  of  unity  in  the  broadside  direction,  0  =  0. 

The  basic  features  of  all  the  minimization  solution  patterns  are  a  placing  of  all 
available  pattern  nulls  in  the  specified  sidelobe  region  and  a  corresponding  degra¬ 
dation  of  the  remainder  of  the  pattern.  For  a  given  number  of  elements,  the 
average  sidelobe  power  in  the  specified  sidelobe  region  increases  with  the  size  of 
the  sidelobe  region  as  the  nulls  become  spaced  more  widely  apart.  For  a  given 
sidelobe  region  the  average  sidelobe  power  in  the  specified  sidelobe  region  decreases 
as  the  number  of  elements  increases. 

We  conclude  that  if  all  except  one  degrees  of  freedom  of  the  element  weights 
are  utilized  to  minimize  pattern  power  in  a  specified  sidelobe  sector,  with  the 
remaining  one  degree  of  freedom  used  to  satisfy  a  mainbeam  constraint  condition, 
then  a  pattern  can  be  formed  with  low  average  sidelobe  power  in  a  broad  sector. 

The  larger  the  number  of  array  elements,  the  lower  will  be  the  average  power  in 
the  specified  sidelobe  sector.  The  price  paid,  however,  for  this  sidelobe  sector 
power  reduction  is  a  marked  degradation  of  the  remainder  of  the  pattern. 
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Figure  5a.  Seven  Element  Pattern  Uniform  Weighting 
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Figure  5b.  Nine  Element  Pattern  Uniform  Weighting 
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